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ABSTRACT 

We present a large, comprehensive survey of rovibrational CO line emission at 4.7 fxin from 69 
protoplanetary disks, obtained with CRIRES on the ESO Very Large Telescope at the highest available 
spectral resolving power (i?=95,000, Aw=3.2kms~^). The CO fundamental band (Aw=l) is a well- 
known tracer of warm gas in the inner, planet-forming regions of gas-rich disks around young stars, 
with the lines formed in the super-heated surfaces of the disks at radii of 0.1-10 AU. Consistent with 
earlier studies, the presence of 100-1000 K CO is found to be ubiquitous around young stars which 
still retain disks. Our high spectral resolution data provide new insight into the kinematics of the 
inner disk gas. The observed line profiles are complex and reveal several different components. Pure 
double-peaked Keplerian profiles are surprisingly uncommon in our sample, beyond the frequency 
expected based on disk inclination. The majority of the profiles are consistent with emission from a 
disk plus a slow (few km s^^) molecular disk wind. This is evidenced by analysis of different classes 
as well as an overall tendency for line profiles to have excess emission on their blue side. The data 
support the notion that thermal molecular winds are common for young disks. Thanks to the high 
spectral resolution, narrow absorption lines and weak emission lines from isotopologues and from 
vibrationally excited levels are readily detected. In general, ^"^CO lines trace cooler gas than the bulk 
^^CO emission and may arise from further out in the disk, as indicated by narrower line profiles. A high 
fraction of the sources show vibrationally excited emission (^50%) which is correlated with accretion 
luminosity, consistent with ultra-violet (UV) fluorescent excitation. Disks around early-type Herbig 
AeBe stars have narrower line profiles, on average, than their lower-mass late-type counterparts, due 
to their increased luminosity. Evolutionary changes in CO are also seen. Removal of the protostellar 
envelope between class I and II results in the disappearance of the strong absorption lines and CO 
ice feature characteristic of class I spectra. However, CO emission from class I and II objects are 
similar in detection frequency, excitation and line shape, indicating that inner disk characteristics are 
established early. 

Subject headings: protoplanetary disks — stars: pre-main-sequence — stars:protostars — 
stars:formation — infrared: general 
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1. INTRODUCTION 

The inner regions of gas-rich protoplanetary disks {R < 
10 AU) arc thought to be the birth places of most gi- 
ant planets (Armitage 2010; Kley & Nelson 2012). The 
chemical and physical processes sculpting these environ- 
ments in the first few million years of the life of a star 
are critical for determining many properties of mature 
planetary systems. These include the formation of rocky 
"oligarchs" - the building blocks of terrestrial planets 
(Nagasawa et al. 2007) - and comets and water-rich as- 
teroids important for the delivery of water and organics 
to planetary surfaces (Raymond et al. 2004). The in- 
teraction between the gas-rich inner disk and protoplan- 
ets has the power to re-arrange the orbital structure of 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 
MS 78, Cambridge, MA 02138; joannabrown@cfa.harvard.edu 

^ Max-Planck-Institut fiir extraterrestrische Physik, Postfach 
1312, 85741 Garching, Germany 

^ Space Telescope Science Institute, 3700 San Martin Drive, Bal- 
timore, MD 21218 

** Leiden Observatory, Leiden University, P.O. Box 9513, NL- 
2300 RA Leiden, The Netherlands 

"" The Kavli Institute for Astronomy and Astrophysics, Peking 
University, Yi He Yuan Lu 5, Hal Dian Qu, Beijing 100871, P. R. 
China 

® Division of Geological & Planetary Sciences, California Insti- 
tute of Technology, Pasadena, CA 91125 

'' ESO, Alonso de Cordova 3107, Casilla 19001, Vitacura, Chile 



the entire planetary system (Armitage 2011), allowing 
for significant modification of radial chemical abundance 
structures; for instance, the chemical boundary defined 
by the snow-line may not always be predictive for the 
compositions of planets in mature systems. 

A key diagnostic of the structure of planet-forming re- 
gions is the fundamental {Av — 1) rovibrational band of 
CO at 4.7 /im. It is particularly sensitive to gas tempera- 
tures of 100-1000 K, corresponding to radii of 0.1-10 AU 
in typical protoplanetary disks around Solar-mass pre- 
main sequence stars. Because of the high opacity of dust 
at 5 /um and the low temperatures of disk midplanes be- 
yond ^ 1 AU, the CO fundamental band typically traces 
the disk at high altitude, specifically the so-called "warm 
molecular layer" (Aikawa et al. 2002; Gorti & HoUen- 
bach 2008; Woitke et al. 2009). In comparison with the 
total disk surface densities in the planet-forming region 
of S = 10 — 1000 gcm^^, the CO fundamental band 
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traces roughly A^h 

to 10~^ — 10~^ gcm"'^, where A^h is the total column of 

hydrogen nuclei. 

While the warm molecular layer represents a small frac- 
tion of the total vertical disk column, its structure is inti- 
mately linked to key properties of the bulk gas in proto- 
planetary disks, relevant to their evolution and ability to 
form planets. Through high temperatures and interac- 
tions with ionizing stellar radiation, the molecular layer 
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acts as a chemical factory, producing water and complex 
organics (Markwick et al. 2002; Glassgold et al. 2009; 
Woods & Willacy 2009; Walsh et al. 2012). It also forms 
a boundary between the deep, neutral and inactive disk 
midplane and the uppermost ionized layers. As such, the 
physical and kinematical structure of the molecular layer 
traces thermal and photo-evaporative flows from the disk 
surface, controlling outward radial mixing and mass loss 
(HoUenbach et al. 1994; Owen et al. 2010). It is in the 
molecular layer that the stellar magnetic field can cou- 
ple to the disk and drive turbulence and accretion flows 
(Gammie 1996; Perez-Becker & Chiang 2011). Finally, 
the molecular layer responds to dynamical perturbations 
caused by the presence of giant protoplanets, leading to 
potentially observable effects (Regaly et al. 2010). 

Rovibrational CO emission is present throughout a 
wide range of protoplanetary disks from still embedded 
protostars (Pontoppidan et al. 2003) to transitional disks 
with inner dust holes (Salyk et al. 2009), and from low 
mass T Tauri stars (Najita et al. 2003) to higher mass 
Herbig stars (Brittain et al. 2007). The CO lines are 
formed by a combination of coUisional excitation, in- 
frared (vibrational) pumping and UV (electronic) fluo- 
rescence. The relative importance of these excitation 
processes depends on location within the disk and the 
shape and strength of the radiation field from the cen- 
tral star (Blake & Boogert 2004). 

The advent of the CRyogenic high-resolution InfraRed 
Echelle Spectrograph (CRIRES) instrument on the Very 
Large Telescope (VLT) has opened the possibility to 
observe the CO fundamental bands at higher spectral 
resolving power (i? ^95000) and higher spatial resolu- 
tion than before, thus providing new insight into this 
critical planet-forming region of the disk. We present 
here the results of a large VLT-CRIRES programme 
of 69 disks around low- and intermediate-mass stars 
and 22 embedded young stellar objects (Pontoppidan 
et al. 2011b). Previous papers have used subsets of our 
CRIRES database to address a variety of questions. Bast 
et al. (2011) investigate the presence and origin of a 
class of broad single-peaked CO rovibrational line profiles 
indicating non-Keplerian gas motions. Using spectro- 
astrometry, Pontoppidan et al. (2008, 2011a) constrain 
the structure of the gas emission and velocity fields on 
milli-arcsecond scales in a smaller sample of disks. Some 
disks show CO emission consistent with simple Keplcrian 
models, whereas other disks, especially those in the Bast 
et al. (2011) sample, show a spectro-astrometry pattern 
consistent with a slow molecular disk wind. For tran- 
sitional disks with a large inner dust hole or gap, Pon- 
toppidan et al. (2008) pinpointed the CO rovibrational 
emission to originate from inside the dust gap whereas 
Brown et al. (2012) resolved the CO emission near the 
outer wall of the hole around the Herbig star Oph IRS 
48. Herczeg et al. (2011) characterized the progenitors 
to protoplanetary disks, when the sources are still em- 
bedded in protostcUar envelopes. Thi et al. (2010) found 
evidence for episodic outflow activity with winds up to 
100 km s^^ in broad blue-shifted lines toward one ob- 
ject, whereas Herczeg et al. (2011) found the same phe- 
nomenon for a handful of other embedded sources. 

Finally, the high spectral resolution of CRIRES also 
boosts the line to continuum ratio and thus allows the 
detection of weak lines from minor species. Smith et al. 



(2009) accurately measured the isotopologue ratios of the 
four major CO species - ^^CO, ^^CO, C^^O and C"0 - 
in the circumstellar environment of two young stars. The 
aim was to search for mass-independent oxygen isotope 
fractionation relative to the interstellar medium in order 
to understand the Solar System oxygen anomaly. Man- 
dell et al. (2012) searched for near-infrared lines of small 
organic molecules (HCN, C2H2, CH4) in a few sources. 
While the these studies discuss the detailed properties of 
specific sub-samples, this paper provides an overview of 
the basic properties of CO rovibrational emission from 
protoplanetary disks around young Solar-type stars at 
the highest available spectral resolution. 

Using the entire CRIRES data set, we aim to address 
the following questions related to the structure and evo- 
lution of planet-forming regions of protoplanetary disks: 
1) Where is the CO gas located? What fraction of sources 
show emission from a radially flowing surface such as 
a slow disk wind? 2) What are typical temperatures 
probed by CO and its isotopologues and how do these 
values compare with current thermo-chemical models of 
the inner disk? How are the CO fundamental bands ex- 
cited, and what consequences does the excitation mecha- 
nism have for the lines as a tracer of inner disk surfaces? 
3) Do CO line profiles depend on stellar spectral type, 
and is this a tracer of the prevalence of UV fluorescent 
excitation over thermal excitation? 4) What evolution- 
ary changes occur in the gas? 

2. OBSERVATIONS AND SAMPLE CHARACTERISTICS 
2.1. The CRIRES survey 

In this paper, we analyze high resolution spectra of the 
4.7 micron CO v=l-0 fundamental emission band using 
the CRIRES on the VLT of the European Southern Ob- 
servatory (ESO) (Kaeufl et al. 2004). CRIRES operates 
at high resolution (R — 95, 000, Av = 3.2 kms"'^) using a 
0'.'2 slit. It is fed by a Multi- Application Curvature Adap- 
tive Optics (MACAO) system, which allows correction of 
atmospheric turbulence and can provide diffraction lim- 
ited images at the focal plane, therefore improving the 
overall instrument sensitivity. 

The sample consists of 91 young stars - 69 protoplan- 
etary disks around young stars and 22 embedded pro- 
tostars that are still surrounded by massive remnant en- 
velopes. The data were obtained as part of an ESO Large 
Programme^ to study infrared molecular emission from 
solar-type protostars and protoplanetary disks (Pontop- 
pidan et al. 2011b)^. We discuss here primarily the disk 
sample, but compare our results with those of the em- 
bedded protostars presented in Herczeg et al. (2011) to 
search for evolutionary trends between embedded and 
classical protoplanetary disks. 

2.2. Sample Selection 

The disk sample spans a range of physical properties 
including spectral type, stellar mass and luminosity, and 
inclination. Such a large sample of CO emission profiles 
provides an opportunity for a broad examination of the 

* This work is based on observations collected at the European 
Southern Observatory Very Large Telescope under programme ID 
179.C-0151. 

^' Fully processed spectra are available on 
http; //www. stsci . edu/~poiitoppi 
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gas distribution in a wide variety of circumstellar envi- 
ronments. The sample was chosen to include protoplan- 
etary disks around solar-mass (0.5 — 2.0 M©) young stars 
(see Figure 1). A range of evolutionary states were sam- 
pled from massive gas-rich disks that are optically thick 
to visible and UV radiation to transition disks with in- 
ner dust holes or gaps. The sources are located in the 
nearby star forming regions visible from Paranal Obser- 
vatory, including Taurus, Ophiuchus, Serpens, Corona 
Australis and Chamaeleon. All these regions have un- 
dergone recent active star formation and have ages of 
approximately 1-5 Myr (Greene & Meyer 1995; Armitage 
et al. 2003; Oliveira et al. 2009). Protoplanetary disks 
were selected based on brightness (> 100 mJy at 4.7 /im), 
as well as the existence of prior datascts characterizing 
the structure of the disks and supporting the presence 
of significant gas reservoirs. The prior datasets defin- 
ing the sample include Spitzer spectroscopy from the 
cores to disks (c2d) Spitzer Legacy survey (Evans et al. 
2003), the Keck-NIRSPEC 3-5 ^im protoplanetary disk 
survey (e.g. Blake & Boogert 2004; Salyk et al. 2011) and 
the VLT-ISAAC protostellar survey (van Dishoeck et al. 
2003) . We limited the survey to disks around Solar-type 
stars (0.1 M© < M* < 2.0 M©), excluding most Herbig 
AeBe stars - the focus of complementary CRIRES sur- 
veys (e.g., van der Plas et al. 2009). 

2.3. Observing procedure and data reduction 

The high spectral resolution and high dynamic range 
of CRIRES spectra, a factor of four improvement in spec- 
tral resolution over most previous CO surveys, fully re- 
solves the individual line profiles so that velocity infor- 
mation can be used to locate the gas within the disk. The 
resolution is particularly needed in the cases of disks with 
low inclinations where the lines are intrinsically narrow 
and in the cases where multiple components, particu- 
larly absorption features, contribute to the line profile. 
Because of the adaptive optics system, line emission can 
be spatially resolved down to angular scales of ^ O'.'l. 

The observations were taken between April 2007 and 
March 2010. Table 1 lists the targets observed and the 
dates and wavelength settings of the observations. Mul- 
tiple spectral settings were taken to cover a range of rota- 
tional J lines. Each spectral setting is observed with four 
different detectors, leading to four discontinuous wave- 
length regions. The wavelength listed is the center of the 
third detector and the range is —70 to -1-39 nm on either 
side of the listed wavelength, with ~6 nm gaps between 
the detectors. The spectra cover the CO P branch lines 
with AJ — —I and a few low energy lines from the R 
branch (AJ — 1). ^^CO v=l~0 lines are the most promi- 
nent. Gas phase lines from higher vibrational states and 
rarer CO isotopologues such as ^^CO, C^^O and C^^O 
are also included in the wavelength settings. 

The spectra were obtained using an ABBA 10" nod- 
ding pattern permitting a first order correction of the in- 
frared background by pair subtraction. Integration times 
were generally 8-16 minutes per spectral setting. The 
reduction of the spectra is described in detail in Pon- 
toppidan et al. (2008, 2011a). Wavelength calibration 
used the prevalent telluric features in the standard star 
spectra. Atmospheric features were removed by divid- 
ing the source spectrum by the standard star spectrum. 
Standard stars were observed close in time and elevation 
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Fig. 1. — A histogram of the spectral types of the stars in our 
disk sample. The sample is dominated by K type stars. When a 
range of spectral types were reported in the literature the mean 
value was taken. 

to minimize atmospheric differences. Remaining small 
discrepancies from airmass differences and subpixel vari- 
ations in wavelength solution were corrected by scaling or 
shifting the standard spectrum to minimize telluric resid- 
uals. Remaining strong telluric features were blanked 
from the analyzed spectra. The observations were pref- 
erentially scheduled for periods when the Earth's veloc- 
ity around the sun created large shifts of the source lines 
relative to the telluric features. When possible, com- 
plete wavelength coverage was obtained using the Earth's 
orbital motion between separated observation dates to 
Dopplcr shift the telluric features relative to the source 
lines. The spectra from two or more dates were then 
combined, barring large differences in line profiles. 

The accuracy of the absolute fiux calibration is ^ 30%, 
estimated by comparing raw counts of consecutive nod 
pairs. This flux calibration is limited by differences in 
Strehl ratio between the target disk and the standard 
star, as well as by pointing jitter. We therefore scaled 
the spectra to known M band fluxes (see Table CI), with 
many from the Spitzer IRAC 2 band. Note, however, 
that some protoplanetary disks are known to be variable 
on short time scales at the 10-60% level (Flaherty et al. 
2012) and discrepancies could exist based on the differ- 
ent aperture sizes. In cases where no Spitzer photometry 
was available a linear relationship between the maximum 
counts received and the flux was determined based on the 
sources with known fluxes. Determining the flux using 
this method led to a median error in derived fluxes of 
~20% for sources with previously known fluxes. How- 
ever, errors were as large as a factor of 2 in exceptional 
cases. The determination was more accurate with larger 
numbers of exposures. 

2.4. Detection rates 

CO gas phase lines are common in our sample, both 
in emission and in absorption, as shown graphically in 
Fig. 2. Lines are commonly a combination of emission 
and absorption resulting from two distinct gas compo- 
nents such as disk emission with foreground absorption. 
Overall, emission features are seen from 53/69 (77%) of 
the disks. Absorption is seen in 44/69 (64 %) of the 
spectra. CO v=2-l lines are seen from 36/69 (52%) of 
the sample, always in emission. ^"^CO emission lines are 
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Fig. 2. — Percentage of disk sources with different lines detected 
in emission (red, horizontally hashed) and absorption (blue, diag- 
onally hashed). In general, rarer isotopologues and higher energy 
lines are less frequently detected. The CO isotopologues are more 
commonly seen in absorption while the vibrationally excited tran- 
sitions are seen exclusively in emission. Table 7 lists common lines 
detected for individual sources. 

seen in 27 (39%) of the sources. However, ^'^CO is more 
commonly seen in absorption (46%), often at cold tem- 
peratures indicative of foreground material. C^^O is seen 
in emission from 4 sources (6%) and in absorption from 
7 (10%). 

2.5. Multiplicity 

In cases where a binary companion was known, we ori- 
ented the slit to observe both stars. There are 8 disk 
systems where both components are well separated and 
bright enough for analysis: AS 205, DoAr 24 E, S CrA, 
SR 24, SR 9, Sz 68, T Tau, VV CrA. Six additional disk 
systems had detectable companions but the M-band sec- 
ondary was too faint for any further analysis: Haro 1- 
4, HD 144432S, HD 144965, SR 9, SX Cha, VW Cha. 
In some cases, the two components are blended in our 
data. HBC 680 consists of a pair of similiar luminosity 
stars with a separation of 0'.'22 (Kohler et al. 2008). The 
northern component of SR 24 is a blend of the B and C 
components with a separation on (/.'OBI (Correia et al. 
2006). 

The spectra from the two components of a binary can 
look very different. AS 205, DoAr 24 E, T Tau and VV 
CrA all have one component seen in emission and the 
other in absorption. In the case of T Tau, the northern 
component is viewed face on while the southern compo- 
nent is seen through its disk (Walter et al. 2003). In 
other cases, one star may be seen through the disk of the 
other (Smith et al. 2009). 



2.6. Extended 
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Good AO correction allows spatially extended sources 
to be seen directly in the CRIRES sample without using 
spectroastrometry. Brown et al. (2012) present a very ex- 
tended ring in IRS 48, whereas Herczeg et al. (2011) ex- 
amine extended emission in two embedded sources. Ex- 
tended emission is not common in our sample, however. 
Only four additional sources have directly imaged ex- 
tended emission: EC 82, LLN 19, R CrA and T CrA 
(see Appendix). Limits for the remaining point-like ob- 
jects are generally around 3 to 4 AU (see also Bast et al. 
2011). 
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Fig. 3. — Gallery of the 8 morphological line profile types dis- 
cussed in Section 3. The profiles are calculated by averaging ob- 
served ^^CO t) = 1 — lines between R(10) and P(32), avoiding 
blended lines. 



2.7. Temporal variability 

Young stars are variable on relatively short timescales 
in a variety of diagnostics from optical fluxes, accre- 
tion line diagnostics and even at longer mid-IR fluxes 
(Flaherty et al. 2012; Muzerolle et al. 2009). Two or 
more epochs in the same setting were obtained for 11 
sources. We approach the search for variability in two 
ways: firstly, by simply comparing the line fluxes from 
the two observations and searching for significant sys- 
tematic differences (taking into account our large error 
bars), and secondly, by closely examining the line profiles 
for differences in shape. 

Overall, we find on timescales of up to two years a 
perhaps surprising lack of variability in most of the line 
profiles. A clear exception is EX Lup, which was un- 
dergoing a large outburst during the period of our pro- 
gramme (see Goto et al. 2011 for further details). RU 
Lup shows a systematic decrease in line flux by 60% be- 
tween April 2007 and April 2008. The decrease could 
be due to either an increase in the continuum level or 
a decrease in the CO flux. Another example is a disap- 
pearance of the VV CrA outflow absorption from 2007 
April to 2007 August. Further details on variability will 
be available in Pontoppidan et al. (2013, in prep). 

3. CLASSIFICATION OF CO LINE PROFILES 

Figure 3 shows representative examples of the different 
CO line profiles of the disk sources. The profiles are often 
shaped by both emission and absorption components. In 
order to estimate the relative importance of the various 
physical components contributing to the observed CO 
lines, we divide the sources into groups based on the 
line profiles. The classification is summarized in Table 
2. We discuss physical interpretations of these different 
categories in Section 5. Figure 4 shows a close up of the 
P(8) line at 4.7359 /im for all sources. In most cases we 
do not find any significant dependence of the emission 
line shape with rotational quantum number within the 
observed range (typically R(l)-P(32)). S CrA B and DC 
CrA are exceptions with a colder narrow component and 
a broader warm component in the ^^CO profile (see also 
Bast et al. 2011). 
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Fig. 4. — A gallery of the P(8) ^^CO line at 4.7359 ^m with the ^^CO R(3) line seen in several spectra at 4.7383 ^ni where available. 
Sources are ordered by similarity of line profile, roughly following the classification scheme in Table 2. RW Aur, SO 411 and WX Cha were 
not covered at these wavelengths and are missing from the plot. 
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3.1. Emission profiles. 

Sources with clean emission lines can be classified 
based on line width (Fig. 5), with four distinct categories: 
(1) narrow, (2) broad single peaked, (3) broad, and (4) 
double peaked. The narrow category sources have line 
widths at base of less that 35 km s~^. Disks in this 
category are known to be face-on and/or transitional in 
nature (TW Hya, LkHa 330, HD 135344 B, IRS 48), with 
CO likely arising from further out in the disk (Salyk et al. 
2011). Lines which are slightly broader (FWHM ~ 13- 
40 km s~^) are generally part of the single-peaked line 
profile sample of Bast et al. (2011). These lines are char- 
acterized by a single central peak with broad wings and 
have large line-to-continuum values. The broad category 
objects retain a single peak but the contrast between the 
line width at the top and bottom (the p value of Bast 
et al. 2011) decreases. The last group shows a double 
peaked structure. 

Many of the CO emission lines in our sample are 
slightly asymmetric. To evaluate the magnitude of any 
asymmetries, the difference between the integrated flux 
in the blue and red sides of the lines was calculated for 
the stacked line profiles (see Figure 6). Where possible, 
the line center was determined from stellar radial veloc- 
ity measurements. However, such measurements do not 
exist for many of the stars in the sample. In these cases, 
the line center was estimated from fitting and checked 
for large deviations from the known molecular cloud ve- 
locities. The lines show an overall tendency for excess 
emission on the blue side, regardless of line profile cate- 
gory. This trend affects lines of all widths indicating that 
it is unlikely to be a geometric effect. Telluric features 
introduce some uncertainty into the integrated fluxes but 
no systematic correlation with the telluric feature place- 
ment is seen. The most dramatic cases of this shift are 
those in the category which show strong emission only 
on the blue side of the line while absorption dominates 
the red side of the line. However, for the majority of the 
sample, the magnitude of the effect is generally small, 
^10% compared to the total ffux. 

3.2. Absorption profiles. 

Absorption lines are common throughout the sample, 
with a detection rate of 67% (46 sources). Line widths 
range from the resolution limit of ~3 km s~^ to tens of 
km s~^. The energy levels with detectable absorption 
also vary from only the lowest J levels to throughout the 
observed range, reflecting underlying temperature differ- 
ences. In general, hotter gas produces broader lines, at a 
level much stronger than expected from thermal broad- 
ening (Fig. 7). The lines occur both with and without 
emission. We divide the lines based on profile shape and 
temperature into four categories: (5) unresolved absorp- 
tion, (6) broad central absorption, (7) blue absorption, 
and (8) absorption with blue emission. The narrow ab- 
sorption lines in category 5 are the most commonly seen 
absorption lines in our disk sample at 59% of all sources 
that show absorption (27/46). For 21 of these sources, 
listed as FA (for foreground absorption) in Table 7, the 
lines are strongest in the lowest J levels and disappear 
for higher energy transitions (J >8). The remaining 6 
sources in this category have hotter CO absorption and 
are listed in Table 2. Warm central absorption (cate- 



gory 6) is less common and almost always occurs with 
resolved line profiles. The > 100 K absorption is gener- 
ally symmetric and centered at the stellar velocity (Fig- 
ure 8) . These warm absorption lines usually have FWHM 
of ^ 10 km s^^ and occur together with broader emission 
lines. Absorption feature velocities generally agree well 
with stellar radial velocities (see Table 3). However, CW 
Tau and SR 24 A and B/C show a shift in the absorption 
line centers redwards for higher energy lines. 

Six sources belong to category 7 with clearly resolved 
blue wings in the absorption lines (Fig. 9). These absorp- 
tion features are usually seen throughout the covered J 
range. The maximum velocities are less than 50 km s~^, 
except for IRS 46 where the gas reaches ^60 km s~^. 
The final category (category 8) have lines that are char- 
acterized by emission on the blue side and absorption on 
the red (Figure 8). 

3.3. Absence of CO lines. 

Only six of our sources show no indication of CO lines 
in either emission or absorption. If CO is present, the 
lines are either weak or narrow enough to be lost in the 
telluric features (<20 km/s). In all cases, the continuum 
is detected (generally S/N of ^^50), indicating the pres- 
ence of circumstellar dust in late-type stars. The CO 
non-detection from UX Tau A is particularly surprising 
as Salyk et al. (2009) show a clear detection within this 
band with NIRSPEC data. However, the telluric feature 
is poorly placed covering -5 to 20 km s~^, which includes 
much of the line center. Also, the lack of CO from the IM 
Lup disk, which has been detected and imaged in CO mil- 
limeter emission out to several hundred AU is puzzling 
(van Kempen et al. 2007; Panic et al. 2009). HD 176386 
is a B9 star (Torres et al. 2006) while HD 144965 A has 
a spectral type of B3 so both disks may simply have no 
CO close to the star due to high photodissociation rates. 
In general, very few sources have no CO visible in the 
spectrum. 

4. CO EXCITATION 
4.1. Rotational temperatures 

The high resolving power of CRIRES allows for the 
detection of low-contrast emission lines, leading to an 
enhanced potential for measuring the rotational temper- 
atures and accurate column densities of the weak opti- 
cally thin isotopologues. In this section, we take advan- 
tage of this to compare the rotational temperatures of 
i^CO, "CO u = 1 - and ^^CO u = 2 - 1 transitions. 

A linear fit to the continuum was determined for each 
line, using uncontaminatcd continuum within 300 km 
s^^ of the line center and 50 km s~^ greater than the 
line width. For accurate measurements of integrated line 
fiuxes, we constructed template line profiles by stacking 
isolated lines free of strong telluric absorption. Absorp- 
tion components were removed using a Gaussian profile. 
The integrated line fluxes were determined by scaling 
the template to each line (or a superposition of tem- 
plates in case of line blends) using Levenberg-Marquardt 
least-squares minimization. This method proved advan- 
tageous to direct integration of the line flux: while direct 
integration works well for strong lines, it is problematic 
for multi-component line shapes, weak lines and lines af- 
fected by telluric absorption. Each line fit was examined 
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Fig. 5. — Normalized ^^CO line profiles of the sources with any emission lines free of absorption, arranged by increasing FWHM. 
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by eye and those that clearly did not recover a meaningful 
line flux due to strong telluric residuals or line blending 
were removed from further analysis. 

The fluxes are put into rotational diagrams to de- 
termine rotational temperatures and column densities. 
Based on an assumed Boltzmann distribution, the ob- 
served flux Fj is such that: 



AttFj 



iVt, 



-E.j/kT„t 



hcj^jQgjAji Q{T^ot) 



(1) 



where VI is the emitting area, Aji is the Einstein A coef- 
ficient of each level i to which the J level can decay, gj 
is the statistical weight of each level (i.e., 2J-I-1), TVtot 
is the column density, Q{T-cot) is the partition function, 
Ej is the energy of the transition, and Trot is the ro- 
tation temperature. For optically thin, isothermal gas 
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Fig. 7. — Comparison of FWHM with rotation temperature for 
absorption components. The resolution limit is marked by the 
dashed line. Unresolved lines generally arise from much colder 
gas than broader absorption lines. The solid line marks the width 
expected from thermal broadening convolved with the instrumental 
resolution. 

in Thermodynamic Equilibrium (TE), rotation diagrams 
produce a straight line with a slope of inverse tempera- 
ture and an intercept proportional to the total mass of 
gas, but are insensitive to emitting area. In the optically 
thick limit, temperature and emitting area are degener- 
ate, while the rotation diagrams are insensitive to the 
total column density. Thus, curved rotation diagrams 
can be a sign of high line optical depths varying with 
J, although multiple temperature components can also 
cause curvature. When both optically thick and thin 
tracers are used, it may be possible to constrain all three 
parameters, iV, Trot and 17, if both isotopologues arise 
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Fig. 8. — Sources with broad absorption lines. Most also have 
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Fig. 9. — Sources with blue absorption wings up to high J, 
attributed to outflows. Most show maximum outflow velocities of 
<30 km s~^ although IRS 46 reaches velocities of ~60 km s^^. 
The VV CrA N data are from 2007 April, by 2007 August the blue 
absorption shoulder is no longer visible. 

from the same region, but, as shown below, this is not 
the case. 
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^^CO is detected in emission from 27 sources and suit- 
able line fits (minumum of good fits to five clean lines) 
could be determined for 12 sources (see Table 4). The 
sources consist almost exclusively of a mixture of transi- 
tional disks and the broad single peaked sources of Bast 
et al. (2011). The lines are relatively weak compared to 
the ^^CO v=l-0 and v=2-l lines with fluxes down by 
factors of ^10 and ^4, respectively. 

The ^"^CO fluxes are distributed along a straight line 
within the errors in the rotational diagrams - consistent 
with single temperature optically thin emission (sec Fig. 
11 for examples). Possible exceptions include AS 205 N 



and DR Tau where some curvature may be apparent in 
the highest and lowest J levels. For the case of DR Tau, 
determination of the ^^CO/C^^O flux ratios indicates op- 
tical depths of the "CO hues of only 0.3±0.2 (Bast et al. 
2011). The emitting radius at which the ^"^CO is opti- 
cally thin is listed in Table 4, assuming a circular slab 
model with uniform temperature and mass. 

Optically thin linear fits uniformly reveal temperatures 
between 200 and 600 K with the majority in the 200-400 
K range (Fig. 10), consistent with the analysis of Bast 
et al. (2011) for a subset of the sources. This is signif- 
icantly lower than has been previously reported for CO 
fundamental band temperatures, generally determined 
from optically thick ^^CO, which may trace different gas. 
Literature values are included in Fig. 10. In the op- 
tically thin case, the range of energies covered should 
have no effect on the linear fit. For the sources with de- 
tectable -^^CO, lines are detected generally up to P/R(15) 
(E„=3656 K/3826 K) and up to R(23) (E„=4722 K) in 
AS 205 N. 

The fact that the ^"^CO lines appear optically thin 
makes it possible to directly correlate the y-intercept and 
total ^■^CO gas mass. Table 4 includes the inferred CO 
column densities assuming an emitting area of 100 AU^ 
(-5 AU radius) and i^CO/i^CO^^eS. The exact area for 
which the ^"^CO emission becomes optically thin depends 
on the specific sources and also on the assumed local line 
profile, taken here to have a total broadening of 2 km 
s^^, the sound speed of H2 at 1000 K, in agreement with 
Herczeg et al. (2011) and Salyk et al. (2011). Larger 
broadening parameters result in more optically thin gas 
while smaller broadening parameters have the opposite 
effect. 

For a few sources, the inferred CO column densi- 
ties are close to 10^° cm'^ (SR 21, LkHa 330) within 
the assumed 5 AU radius, which would imply values of 
Nh > 10^^ cm~^ for a standard conversion factor of 
CO/H2=10~'*. These column desities are very large com- 
pared to the surface column densities of Nh ^ 10^^-10^'^ 
cm~^ that are expected down to the layer where the 5 /im 
continuum becomes optically thick (e.g., Aikawa et al. 
2002; Gorti & HoUcnbach 2008; Woitke et al. 2009). An- 
other way to look at this problem is to derive total gas 
masses from the ^^CO data. Values computed under the 
optically thin assumption (but independent from any as- 
sumed emitting area) are included in Table 4. For some 
sources, the masses are >10~^ M0, implying that a sig- 
nificant fraction of the disk mass would be contained in 
just the surface layers of the inner few AU of the disk, 
under these assumptions. One possible solution to this 
conundrum of the large inferred column densities and 
masses is that the dust grains are settled to the mid- 
plane so that the gas/dust ratio is much larger than the 
standard value, which allows us to look deeper into the 
disk at 5 iim. A correlation between mid-IR SED slopes, 
a tracer of settling, and CO equivalent widths indicates 
that dust settling may increase the observable CO (see 
§5.4). Another possibility is that UV or IR radiative ex- 
citation and resonant scattering contributes to the ^^CO 
lines fluxes (see §4.1.2), in which case the inferred column 
densities and masses are upper limits. 

The results from optically thin flts of the ^'^CO gas, 
detected here with a much larger frequency than in pre- 
vious studies, imply a different location for at least some 
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Fig. 10. — (left) ^'^CO temperatures derived from optically thin fits of disk emission sources compared to temperatures derived from 
fitting the ^■^CO lines, (right) ^^CO temperatures derived from optically thin fits compared to CO temperatures reported in the literature 
(Salyk et al. 2009, 2011). The three literature sources with ^■^CO temperatures in the 200-300 K range are all sources for which ^^CO was 
previously detected. 

vat ions. 

In theory, optically thin ^^CO can constrain the degen- 
eracies in modeling the optically thick ^^CO emission, 
assuming that the emission arises from the same gas. 
However, the low ^"^CO temperatures are incapable of 
explaining the observed ^^CO high J lines in most cases, 
indicating the presence of a warm component seen in the 
i^CO (Fig. 11). We examine HD 135344 B in detail to 
determine whether a two temperature model is capable 
of fitting this behavior (Fig. 12). We use a simple slab 
model which calculates the optical depth in each transi- 
tion (see e.g., Salyk et al. 2009; Brown 2008). The "CO 
data are best fitted with T,ot=250 K, iV(i3C0)= 8 x lO^^ 
cm~^ and an emitting region with a radius of 5 AU to en- 
sure that the ^■^CO emission is optically thin. Using such 
a large radius and multiplying the column density by 65, 
the low ^'^CO temperature is incapable of explaining the 
high J P (20-27) ^^CO lines which are clearly detected 
(formal errors are smaller than the points) (Fig. 11). A 
smaller, hotter (850 K) region of ^^CO emission needs to 
be added to fit the ^^CO data. We conclude that a two 
temperature model (black line) is capable of fitting the 
data as the expected ^'^CO fluxes from the warm ^^CO 
component are below our detection threshholds if close 
to optically thin. 

Previous studies have also had difficulty reconciling 
the ^^CO and ^^CO line fluxes for some sources. Salyk 
et al. (2009) found a similarly low temperature of 250- 
350 K for I^CO in SR 21. This is compatible with the 
^^CO lines due to a lack of detected high J lines for that 
source. However, the required mass of 6x10^^ g of CO 
(gas mass of 3 x 10~^ M©) and emitting area of 78 AU^ 
are large, although a gas location at radii as large as 
7 AU is confirmed by spectro-astrometry (Pontoppidan 
et al. 2011a). Blake & Boogert (2004) were unable to 
explain the high relative fluxes from ^''CO in AB Auri- 
gae and Brittain et al. (2009) found anomalously high 
^■^CO line fluxes from the disk around A-type star HD 
100546, with an apparent ratio of ^^CO/"CO of 4 rather 



Fig. 11. — Rotational diagrams of ^^CO (filled circles) and ^^CO 
(open circles). A linear fit was made to the ^''CO, lower dashed 
line, and was scaled up by an isotope ratio of 65 to compare to 
the expected ^^CO if optically thin. The high J lines are clearly 
significantly underpredicted and require a hotter gas component 
to fit. The lower J lines are likely underpredicted due to optical 
depth effects. 

of the CO than indicated from studies of ^^CO. The low 
temperatures and large emitting areas suggest that much 
of the ^'^CO arises from either larger disk radii or deeper 
into the disk. 
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The ^^CO lines throughout our sample show curvature 
in the rotation diagram which has been attributed to op- 
tical depth effects (Blake & Boogert 2004). They look 
qualitatively similar to previously published rotation di- 
agrams from disks (Salyk et al. 2009; Najita et al. 2003). 
In cases where the sources had been previously observed 
at lower spectral resolution (e.g. NIRSPEC, Salyk et al. 
2009), benchmark tests ensured that our CRIRES rota- 
tion diagrams quantitatively agree with previous obser- 
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Fig. 12. — The rotational diagram for ^-^CO (filled circles) and 
^^CO (open circles) from HD 135344 B. The blue and purple lines 
are isothermal slab models based on the rotational temperature 
and mass determined from the optically thin fits to the ^^CO (250 
K, 8 10^^ cm"-^) for two emitting radii. The yellow line is a higher 
temperature model (850 K, A''(^^CO)=10^^-® cm~^ column den- 
sity) to fit the high J ^■^CO and the red line is this model combined 
with the 5 AU radius ^^CO model. 

than the interstellar ratio of ^65. In both cases, non- 
thermal processes were invoked to explain the discrep- 
ancy. Blake & Boogert (2004) find that resonant scat- 
tering of a large fraction of the IR continuum photons can 
reproduce the AB Aur fluxes. Brittain et al. (2009), on 
the other hand, suggest that UV pumping is the primary 
excitation mechanism and the overabundance of ^'^CO is 
the result of additional sources of opacity in the disk. 

IR resonant scattering (IR pumping) of ^^CO primarily 
enhances the low energy J lines, reflecting the thermal 
distribution of the lower vibrational level population in 
the excited state. Optical depth effects may reduce the 
strength of this component in the ^^CO lines. Excitation 
via IR pumping is expected to be particularly effective 
for disk wind sources and transitional disks, explaining 
why we see ^^CO mostly from these sources. Winds lift 
material above the disk where the CO molecules are more 
exposed to radiation. In transitional disks, the lower 5 
yum continuum optical depth may result in less shielding 
from the dust and thus enhanced absorption. Resonant 
scattering is primarily dependant on the received flux 
and would imply column densities several orders of mag- 
nitude smaller to reproduce our observed line fluxes. 

4.2. Vibrational excitation 

Vibrational lines from CO v=2-l, v=3-2 and v=4-3 
are seen from a subset of sources. CO v=2-l lines are 
the strongest and thus most commonly detected. These 
higher energy lines are likely pumped by UV flux as the 
upper level temperatures of ^7000 K are too high for 
coUisional excitation by the 500-1000 K gas seen in the 
i^CO lines. 

A vibrational flux ratio was calculated for each source 
based on the v=2-l flux divided by the v=l-0 flux for 
the equivalent rotational line (Table 5. UV pumping does 
not induce large changes in population distribution in the 
different rotational states since the selection rule only al- 



lows quantum number changes of AJ ±1. Thus, the 
shape of the v=2-l rotational diagram is similar to the 
v=l-0 lines and the offset between the two reflects the 
difference in upper state population. Measuring the vi- 
brational flux ratios of equivalent J rotation lines thus 
traces the population difference. Typically five flux ratios 
were found and the median was taken for the final value 
to remove outliers. The vibrational flux ratios range from 
0.05 to 0.5 and upper limits were calculated for high S/N 
spectra where the CO v=2-l lines were not detected. 

The observed vibrational flux ratio generally increases 
with accretion luminosity (Fig. 13). The accretion lu- 
minosity is determined based on simultaneous measure- 
ments of the Pf/3 line which is covered serendipitously in 
the CRIRES spectrum using the conversion to accretion 
luminosity as given in Salyk et al. (2013). IRS 48 is an 
outlier in the plot, probably due to its higher spectral 
type of AG resulting in a photospheric UV flux rather 
than accretion luminosity pumping the v=2-l lines and 
the 30 AU gas location (Brown et al. 2012) reducing col- 
lisional excitation, and is therefore excluded from the 
linear fit. 

We modeled the effects of a stellar UV field on the 
CO emission in a UV excitation model similar to Brit- 
tain et al. (2007). UV fluorescent pumping of the CO 
vibrational lines in the fundamental band involves exci- 
tation by UV photons from the X^T,'^ ground electronic 
state to the A^H excited electronic state. The decay 
from the excited electronic state to the ground state can 
increase the vibrational energy of the molecule by pop- 
ulating higher vibrational levels in the ground electronic 
state than would be expected from purely thermal exci- 
tation. We assume that the system is in a steady state 
and balance pumping into the ^^11 state from the inci- 
dent UV radiation field with spontaneous emission out of 
the A^Ii state. Spontaneous emission within the X^I]+ 
ground state vibrational levels is also included. We ex- 
panded the vibrational states included in the model to 35 
X^I]+ ground levels and 25 ^^11 levels; rotational levels 
are not explicitly taken into account. For this purpose, 
we computed the oscillator strengths and Einstein-A co- 
efficients between individual A-X vibrational levels us- 
ing the Rydberg-Klein-Rees (RKR) program of Le Roy 
(2004). Accurate potential curves based on the spectro- 
scopic data of Le Floch (1992) for the X state and Field 
et al. (1972) for the A state were used together with the 
A-X transition dipole moment of Gilijamse et al. (2007). 

One of the main reasons for the expansion to higher 
vibrational states was to probe the effects of Lyman a 
emission (1216 A), which dominates the far ultraviolet 
emission from classical T Tauri stars (Schindhelm et al. 
2012) and can pump the v=14-0 line of CO (France et al. 
2011). Much of the observed Lyman a is absorbed by the 
ISM before reaching Earth so the flux seen by the disk 
must be reconstructed from fluoresced H2 lines (Herczeg 
et al. 2002). However, even with reconstructed Lya line 
profiles, the effects of Lyman a are only noticeable at 
high vibrational levels but are marginal at the v=2-l 
level that is seen in our sample. Some of this result can 
be explained by the fact that the oscillator strength for 
the 14-0 band is more than three orders of magnitude 
smaller than those for the 1-0, 2-0 and 3-0 bands which 
dominate the UV pumping. The other factor is that the 
downward relaxation in the X band is much faster than 
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Fig. 13. — Flux ratio in vibrationally excited levels derived 
from the 2-1/1—0 flux ratios versus accretion luminosity deter- 
mined from simultaneous Pf/3 emission. The accretion luminosity 
in the UV drives UV fluorescent pumping which is reflected in in- 
creased CO v=2-l line fluxes. IRS 48 is an outlier in the plot likely 
due to its higher spectral type of AO resulting in photospheric UV 
flux rather than accretion luminosity pumping the v=2-l lines, 
and is therefore excluded from the linear flt. 

the UV excitation, spreading the excess flux through all 
the lower energy vibrational levels. 

Fig. 14 presents our model results for the vibrational 
population fraction of the v=2/v=l vibrational levels 
for a pure UV fluorescent pumping model using differ- 
ent blackbody radiation fields. The vibrational popula- 
tion fraction decreases with stellar effective temperature 
due to the smaller amount of UV available to pump the 
higher vibrational levels of the A state. The observed 
vibrational flux ratios are included in Fig. 14 and gen- 
erally show much greater values than would be expected 
from simple blackbodies, especially for the later spectral 
types. This is in agreement with the correlation found 
with accretion luminosity in Fig. 13. For the cases of 
TW Hya and DR Tau, both strongly accreting objects, 
we have also run models using the observed stellar spec- 
tra including the enhanced UV emission (Herczeg et al. 
2002; Yang et al. 2012), resulting in vibrational flux ra- 
tios due to pure UV pumping of ~0.4-0.5. The observed 
flux ratios for these sources are 0.05 and 0.2, respectively, 
indicated with blue dots in Figure 14. The difference be- 
tween models and observations is likely due to excess 
population in the v=l levels from IR pumping or col- 
lisional excition, which results in a lower vibration flux 
ratio. More detailed modeling and observations of higher 
energy vibrational lines are needed to disentangle these 
two excitation processes (e.g. Brittain et al. 2007, Brown 
et al. 2012, Bast et al. 2011). However, the ubiquity and 
strength of v=2-l lines in the sample appear consistent 
with UV fluorescent pumping primarily from accretion 
luminosity. 

5. DISCUSSION 

5.1. Where is the CO gas located? 

The presence of multiple line profile components and 
temperatures suggests that separate, physically unre- 
lated gas reservoirs are contributing to the spectra. Such 



Fig. 14. — Vibrational flux ratio (v=2-l/v=l-0) versus spec- 
tral type compared to the pure fluorescence model for vibrational 
emission (dashed line). Filled black circles are the observed values, 
open circles are model results for blackbodies at 4000, 6000, 8000 
and 10000 K and squares are model results using the observed UV 
spectra of DR Tau and TW Hya while the blue dots are the ob- 
served values. The observed increased UV flux over blackbody is 
clearly essential at late spectral types. 

reservoirs could include 1) superheated gas from the in- 
ner (< 1 AU) disk surface (Najita et al. 2003; Blake & 
Boogert 2004), 2) emission/absorption from an extended 
disk wind or other outflow activity (Pontoppidan et al. 
2011a; Herczeg et al. 2011), 3) absorption from the outer 
disk in an inclined system (Rettig et al. 2006), 4) ab- 
sorption from another disk around a binary companion 
(Smith et al. 2009), 5) absorption from remnant enve- 
lope (Boogert et al. 2002a; Herczeg et al. 2011), and 6) 
absorption from cold foreground molecular gas, entirely 
unrelated to the young star (Boogert et al. 2002b). 

Inner disk: The surface layers of disks have long been 
thought to be one of the primary producers of funda- 
mental CO emission (e.g. Najita et al. 2003). The gas 
in the disk atmosphere has regions with the 100-1000 K 
temperature range needed to thermally excite the v=l 
level of CO and produce the abundant strong CO v=l-0 
lines. The line widths at the base are consistent with 
expected Keplerian rotation in the inner regions of disks 
(Salyk et al. 2011, §B). 

Disks with a pure Keplerian velocity field (no radial 
motions and Vphi oc R^^''^) produce a characteristic line 
profile based on inclination. A double peak arises from 
inclined disks due to the radial decrease in temperature 
and resulting decrease in low-velocity CO emission from 
the cold outer disk causing a central dip. The double 
peaked line profiles (category 4) are good examples of this 
classic Keplerian profile and can be seen in the cases of 
GQ Lup and RNO 90. The central dip disappears at low 
inclinations due to the lack of velocity contrast between 
the inner and outer disk resulting in narrow single peaked 
lines (category 1), such as TW Hya and HD 135344 B. 
In these cases, the two peaks can still be separated using 
spectro-astrometry (Pontoppidan et al. 2008). 

The list of sources with clear classical Keplerian line 
profiles consistent with these models is given in Table 2 
and is surprisingly small. As can be seen in Figure 5, 
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the majority of our line profiles are single peaked: our 
percentage of double peaked sources is only 10% (7/69) 
with 6 additional sources (category 6) which may be 
self-absorbed highly-inclined disks. This lack of double 
peaked profiles has been noted in even the earliest CO 
surveys (Najita et al. 2003). However, the constraints 
on Keplerian models are greater with the improved spec- 
tral and spatial resolution and much large sample. If the 
inclinations of the full CRIRES sample are randomly dis- 
tributed and all line profiles are taken to be Keplerian, 
these low statistics would imply an inclination of >70° 
before a double peaked line profile would be resolved. 
For a 1 Mq star at the CRIRES spectral resolution, this 
would imply some CO emission from ~10 AU or greater, 
which is not seen in the 2-D spectral traces. For less 
massive stars, limits on the outer extent are less extreme 
(e.g. 0.5 M0 corresponds to ^5 AU). However, no trend 
is seen with spectral type and uniformly low masses are 
excluded by the range of measured stellar types. We 
conclude that the warm surface gas in the inner disks 
in the majority of our sources has kinematics that are 
inconsistent with only Keplerian rotation. 

Disk winds: Many of the sources with the strongest line 
to continuum ratios are single peaked with broad wings 
(category 2), including AS 205 A, DR Tau and RU Lup. 
As discussed in Bast et al. (2011), the combination of 
a narrow single peak and broad wings is incompatible 
with pure Keplerian profiles without high temperatures 
in the outer disk. The lack of spatial extent rules out 
this model. Indeed, Pontoppidan et al. (2011a) use spec- 
troastrometry to show that the spatial distribution of the 
emission is far more compact and asymmetric than can 
be explained by such a Keplerian disk. They propose a 
slow moving disk wind as the source of the central emis- 
sion, but with the broad line wings still due to Keplerian 
rotation imparted from the launching region of the disk. 
The line profiles appear symmetric and centered close to 
the source radial velocity (within ^ 5 km s~^, see Table 
3) constraining the physical characteristics of the wind. 
Sources with these profiles show higher vibrational CO 
lines indicating fluorescent excitation by the UV radia- 
tion produced by the accretion shocks. 

The largest category of emission lines is category 3 
consisting of the broad profiles (> 200 km s~^). These 
were not included in the Bast et al. (2011) sample due to 
lower contrast between the line width at 10% compared 
to 90% of peak height and/or lower line to continuum 
ratios. However, based on the prevalence of these single 
peaked lines, they may have the same origin in a wind. 

Outer disk: Category 6 (emission with broad central 
absorption) sources have strong absorption lines above 
P(10) (4.7545 ^J,m, E^ppcr=3330 K) and likely arise from 
the surface layers of the outer disk. These absorption 
lines commonly occur in conjunction with disk emission 
lines although the emission can be difficult to categorize 
due to the absorption in the line center. Examples of 
warm CO absorption occur in CW Tau, SR 24, and T 
Tau S. Temperatures range from 100 to 300 K (§ 4). 
At near edge-on inclinations, the outer disk atmosphere 
could absorb CO as has been suggested in the case of 
T Tau S (Rettig et al. 2006). In the case of binaries, 
it is possible that the absorption can arise from looking 
at the star through the companion's disk as may be the 
case in the VV CrA system where VV CrA N would be 



seen through the VV CrA S disk (Smith et al. 2009). 

Outflows: Molecular outflows and winds moving to- 
wards us in our line of sight result in absorption lines 
with broad blue-shifted wings seen in category 7 sources. 
The gas producing the lines is warm (up to 1000 K) likely 
originating from currently shocked gas rather than en- 
trained outflow gas. The lines are blue shifted due to 
the disk blocking the red-shifted component and its po- 
sition behind the illuminating star. The wings are not 
always smooth indicating clumps in the wind as in the 
case of DC Tau (shown in Figure 3). LLN 19 has episodic 
outflows (Thi et al. 2010). VV CrA N also has outflow 
components which disappear between April and August 
2007. The profiles are similar to those found for the em- 
bedded sources in this CRIRES sample highlighted in 
Herczeg et al. (2011). 

The category 8 sources with blue emission and red ab- 
sorption are probably from a combination of outflow in 
emission and residual envelope in absorption. Both the 
emission and absorption have similar warm temperatures 
of a few hundred K. These lines are the most asymmet- 
ric of the categories. The warm absorption is close to 
the cloud velocity and the emission could then be from 
an outflow. In the case of T CrA, which is known to 
drive a jet (Wang et al. 2004), the emission component 
is spatially extended deflnitively indicating an origin in 
an outflow. 

Foreground molecular clouds: Many of the CRIRES 
sources are towards high extinction regions of nearby 
star forming clouds. CO gas should be present in the 
foreground molecular cloud with low turbulent broad- 
ening and cold temperatures (~10-100 K). We attribute 
the Category 5 narrow unresolved absorption lines which 
appear only at low J levels to this foreground material. 
This combination of cold temperatures and low densities 
is conflrmed by rotational diagrams (see § 4) . These lines 
are very common throughout the sample with a detection 
rate of 38%. The high spectral resolution of CRIRES is 
a strong advantage for complex line profiles because the 
foreground lines can be readily separated from emission 
components and other dips. The binary pairs, such as 
DoAr 24E also show the dumpiness of foreground ma- 
terial with even cold absorption components varying be- 
tween binary pairs on spatial scales of a few 100 AU. 

5.2. Temperature structure of the inner disk region 

Typical temperatures of the ^'^CO lines appear lower 
(<500 K) than has generally been derived from likely 
optically thick ^^CO lines (500-1000 K). ^^CO is close to 
being optically thin and thus may be able to probe deeper 
into the disk's vertical structure than ^^CO. The lower 
temperatures indicate that the ^"^CO lines are dominated 
by emission from either further out or deeper vertically 
in the disk. As discussed in §4.1.1 and 4.1.2, both UV 
and IR pumping may contribute to the line fluxes, but 
these processes should not change the inferred rotational 
temperatures significantly. 

Line widths systematically differ between ^^CO v=l- 
0, v=2-l and ^^CO, indicating that the lines do not arise 
from identical locations in the disk. For each disk, the 
CO v=2-l line profiles are generally broader than v=l-0 
while the ^^CO lines are generally narrower than ^^CO 
v=l-0 (Fig. 15). This trend was noted by Bast et al. 
(2011) in their broad single peaked sources but appears 
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Fig. 15.— "CO FWHM/12CO FWHM (dots) and i^co v=2-l 
FWHM/l^CO v=l-0 FWHM (open circles) versus FWHM of the 
v=l-0 lines. The ^^CO v=2-l line widths range from similar to the 
v=l-0 lines to much broader. The ^'^CO line widths are generally 
down by a factor of ~2. The two sources with similar FWHM 
are SR 21, where all emission lines appear to arise from the same 
region, and S CrA S, which has two distinct components in the 
^^CO with the cooler ^^CO component having the same line shape 
as the "CO. 



more generally throughout the sample. The change in 
average width between the ^^CO and ^"^CO lines corre- 
sponds to a change in radius by a factor of ~4, assum- 
ing Keplerian rotation. If the entire temperature differ- 
ence between the isotopes is assumed to come from a 
radial difference, this would imply a temperature profile 
T ex RP where p is between -0.5 and -1.5. 

We can also compare our temperature and thermal col- 
umn density measurements with current chemical mod- 
els. Many thermo-chemical models of inner disk regions 
exist and they typically have gas kinetic temperatures of 
1000 K or larger in the layer where CO becomes abun- 
dant (e.g., Gorti & Hollenbach 2008; Gorti et al. 2011; 
Glassgold et al. 2009; Woitke et al. 2009; Najita et al. 
2011; Walsh et al. 2012; Bruderer et al. 2012). Ahhough 
there is a large spread in the model results, these tem- 
peratures are broadly consistent with those measured for 
i^CO. 

Woitke et al. (2009), their Figure 14, provide gas tem- 
peratures weighted by the CO abundance as functions 
of CO column density for different disk radii. Their 
model indicates that the low J CO rovibrational lines 
become optically thick at 1 AU for a column density of 
10^^ cm~^, at which point the CO weighted temperature 
is 600 K. Rarer isotopes, such as "'^'^CO, are sensitive to 
larger column densities before becoming optically thick 
and thus probe lower temperatures. Our CRIRES ^"^CO 
measurements follow the predicted relation between col- 
umn density and gas temperature in their Figure 14 but 
lie above the 1 AU contour by a factor of ^3x10^ in col- 
umn density. This may be further indication that the 
^■^CO lines are not just thermally excited and thus the 
column densities are overestimated. 

Gorti et al. (2011), their Figure 8, model the CO rovi- 



FlG. 16. — A plot of the line width at base vs spectral type. 
This figure shows that earlier type stars tend to have narrower 
line widths. Our CRIRES sample (filled circles for late-type stars, 
crossed circles for early type) consists mainly of later type stars. 
Additional Hcrbig stars from Brittain et al. (2007, 2009) and van 
der Plas ct al. (2009) are included as open circles. 

brational lines in the transitional TW Hya disk and find 
that inner disk masses between 10^^ to 10~^ Mq fit the 
rovibrational lines well to within a factor of 2. Most of 
the model emission is dominated by the warm surface 
layers rather than deeper material. This amount of in- 
ner disk gas is similar to those derived from our ^"^CO 
data (Table 4). 

5.3. Effects of stellar mass on CO 

The fundamental band of CO from Herbig AeBe stars 
has been studied extensively (Najita et al. 2003; Blake 
& Boogert 2004; Brittain et al. 2007; van der Plas et al. 
2009) . One of the most obvious differences is the preva- 
lence of high vibration lines in Herbig stars with transi- 
tions up to v=8-7 being seen (Brittain et al. 2009). In 
comparison, the T Tauri stars in our sample show vibra- 
tionally excited emission generally only to v=2-l with a 
significant fraction showing no vibrationally excited lines 
at all. This is probably due to the weaker UV fields emit- 
ted from T Tauri stars as is clearly predicted from models 
(Figure 14). 

A comparison of the widths of the line bases (see Ap- 
pendix) reveals that Herbig stars tend to have narrower 
lines (Figure 16). This indicates that the inner edge of 
the gas is further out in the disk compared to T Tauri 
stars. As our sample consists mainly of later type stars, 
additional Herbig stars from Brittain et al. (2007, 2009) 
and van der Plas et al. (2009) are included. Our few Her- 
big stars show similar line widths indicating that this is 
unlikely to be due to any differences in methodology. A 
Kolmogorov-Smirnov test places the probability at 99% 
that the samples have different underlying distributions. 

For a given temperature, the radial location will be 
further out in a Herbig disk than in a T Tauri disk due 
to the larger luminosities of Herbig stars (Salyk et al. 
2011). As the CO lines are sensitive to only a certain 
range of gas temperatures, the lines would be expected 
to arise from further out in Hcrbig disks. While this 
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Fig. 17. — A comparison of the spectra from the class I source IRS 43 (top) and the class II disk RNO 90 (bottom). The deep ice feature 
at 4.67 /tm is only present in the class I source. RNO 90 has absorption lines only in the lowest J ^^CO lines from cold foreground gas. 
Most of the differences between the two spectra can be attributed to the disappearance of the envelope. 

the outer disk. 

The emission hnc characteristics of the class I sources 
are similar to those of the class II sources (see Herczeg 
et al. 2011 for further details). The class I profiles are 
complicated by absorption from the envelope and out- 
flow but the underlying emission seems similar to the 
class II disk emission with the addition of more common 
blueshifted components. In particular, the class I emis- 
sion lines have similar widths, temperatures and lack of 
spatial extent to the class II sources. 8/16 of class I 
sources show CO v=2-l emission lines (although 3 show 
no ^^CO in emission) which is similar to the 47% of class 
II sources which show CO v=2-l emission. Unlike the 
class II sources, the class I sources show no double peaked 
line profiles at all, but are always single-peaked indica- 
tive of a slow disk wind. Overall, the inner few AU of 
class I and II disks are very similar indicating that disk 
characteristics are established early. 

Full to transitional disks: Transitional disks with inner 
gaps and holes provide an opportunity to study inner disk 
gas at a time when the dust is disappearing. The tran- 
sitional disks generally fall into the narrow line profile 
category, partly coincident with the known close to face- 
on inclinations of many of the disks and partly due to the 
generally narrower line widths of transitional disks found 
by Salyk et al. (2011). The lower dust continuum should 
increase the effectiveness of radiative excitation in the 
inner regions. This appears to be the case with all tran- 
sitional disks showing UV-pumped CO v=2-l emission 
and most showing the rarer ^^CO emission lines which 
may be IR pumped. 

SED slopes provide clues to disk structure and evolu- 
tionary states. We focus on the ratio between 13 and 30 
^im (Fix3o,2m/Fi/i3/im) as defined in Brown et al. (2007). 
Low numbers likely indicate strong settling while ex- 
tremely high numbers mark transitional disks with inner 
gaps and holes (see Figure 19). A trend is seen with the 
lowest ?T.i3_3o sources having the strongest P(8) equiva- 



explains the general trend that we see, it does not explain 
the prevalence of Herbig stars with large (> 10 AU) gas 
holes (IRS 48 in our sample, Brown et al. 2012; see also 
Goto et al. 2006, van der Plas et al. 2009, Brittain et al. 
2009). Some additional effect, e.g. photodissociation of 
CO, is needed to explain these holes. 

5.4. Evolutionary changes in CO 

Complexity in line profiles generally decreases at each 
evolutionary stage. Transition disks thus show the most 
straightforward profiles and embedded class I sources the 
most kinematically complex. 

Class I to Class II: The 18 class I embedded objects 
discussed by Herczeg et al. (2011) show some system- 
atic differences from the 69 disk sources. The first is 
that the CO ice band seen prominently in class I sources 
(Pontoppidan et al. 2003) disappears from the spectra of 
the disks due to the disappearance of the cold envelope 
(Fig. 17, 18). There are only a few exceptions of class II 
disks with ice bands, with the ice most likely arising in 
dense foreground gas. Disk winds/outfiows as refiected in 
blue absorption wings are seen in both classes. However, 
these sources are less common in the disk sample and, 
when seen, the maximum velocities are smaller (gener- 
ally <20 km s~^ as opposed to 50-100 km s^^). 

The number of sources with gas phase absorption lines 
decreases from 100% to ^60% as sources evolve from the 
embedded class I to the class II phase (Fig. 18). A com- 
parison of the absorption line temperatures reveals that 
the class II sources can show very low temperature (<50 
K) components (Fig. 7). However, both class I and class 
II sources can have warm absorption lines. The inner 
regions of class I envelopes probably result in higher av- 
erage temperatures along the line of sight while class II 
sources may have absorption lines arising solely from cold 
foreground material. The few class II sources with hun- 
dred K absorption lines likely have an inclination angle 
such that absorption arises through the surface layers in 
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Fig. 19. — A trend is seen for declining SEDs, as seen in the 
13-30 ^m flux ratios, to have larger CO P(8) equivalent widths. 
The flux ratios are taken from Brown et al. (2007); Furlan et al. 
(2009), and the Spitzer archive. The correlation has a Pearson's 
correlation coefficient of -0.75 indicating a strong correlation. 

lent widths and the highest fT.i3_3o sources having weaker 
P(8) equivalent widths. This correlation is not reflected 
in the fluxes, indicating that the cause lies in the inter- 
play between dust continuum and gas lines. This may be 
an optical depth effect where more CO resides above the 
dust T — I surface in the settled disks leading to strong 
CO lines. 

6. CONCLUSIONS 

CO fundamental lines are commonly seen from the cir- 
cumstellar environments of young stars and provide in- 
formation on molecular gas in the planet-forming region 
in disks that ALMA will not be able to probe. Our ex- 



tensive collection of high spectral resolution CRIRES ob- 
servations reveal complex environments with signatures 
of disks, envelopes, foreground molecular clouds, winds 
and outflows in diverse combinations. 

Young disks display a variety of complex line profiles 
with simple double-peaked Keplerian disk profiles being 
surprisingly rare. Most line profiles have wings consistent 
with Keplerian rotation in the inner disk close to the 
dust sublimation radius. However, the majority of the 
emission lines have excess flux at the line center, which in 
a disk model requires emission from large radii. In most 
cases, there is no evidence of spatial extent in the CO 
lines, ruling out standard Keplerian disk models. This 
trend, plus a tendency for slight excess emission on the 
blue side of the lines, leads us to conclude that slow disks 
winds may be common from young stars and traceable 
through CO lines. 

The high spectral resolution and sensitivity of CRIRES 
allow a more comprehensive overview of emission lines 
from weaker isotopologues. ^■^CO lines suffer less con- 
fusion from optical depth effects than ^^CO lines and 
thus provide better constraints on the gas properties. 
The ^"^CO lines are detected with much greater frequency 
than in previous studies and have lower excitation tem- 
peratures, indicating an origin further out radially or 
deeper within the vertical structure. The ^'^CO line 
widths are generally narrower than those of ^^CO, in- 
dicating dominance by slower velocity gas arising from 
larger radii. The ^'^CO line strengths for the objects with 
the coldest gas imply large column densities/emitting 
areas if only collisional excitation is assumed. IR res- 
onant scattering may contribute to the line fluxes, as 
has been proposed for Herbig Ae/Be stars, so that the 
inferred column densities are upper limits. The de- 
rived temperatures are roughly consistent with current 
thermo/chemical models of the inner disk. 

Non-thermal excitation is clearly reflected in the CO 
lines, hinting at the complex photoprocesses at work in 
the molecular layer. Vibrational emission from thermally 
inaccessible energy levels is detected in about 50% of 
the sources. This emission is correlated with accretion 
luminosity indicating that the UV radiation arising from 
accretion drives the fluorescent pumping, particularly in 
late-type stars where the stellar UV continuum is weak. 

Finally, we examine the large sample for trends be- 
tween different types of sources. The change from class 
1 embedded protostar to class 11 disk produces dramatic 
changes in the CO spectra, as the removal of the proto- 
stellar envelope results in the disappearance of the strong 
absorption lines and ice feature characteristic of class 1 
spectra. However, the emission lines from class 1 and 11 
sources, both arising from close to the star, are similar in 
detection frequency, excitation and line shape, indicat- 
ing that disk characteristics are established early. The 
beginning of disk dispersal seen in transition disks also 
brings changes to the CO spectra. Transition disk CO 
lines are generally Keplerian with a tendency for narrow 
line widths indicating gas at larger radii. Herbig Ae/Be 
stars also typically have narrower lines, but in this case, 
it is likely due to an increase in the dust sublimation ra- 
dius and the disk region with temperatures suitable for 
detecting rovibrational CO. 

We have highlighted general trends from the large sam- 
ple to provide an overview of CO gas within young sys- 
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terns. However, the individual spectra are all unique, 
reflecting physical differences among the stars, disks and 
surrounding environments. We hope that future detailed 
modeling will be able to exploit fully all the information 
contained in these rich CO spectra. 
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APPENDIX 
SOURCES WITH EXTENDED EMISSION 

Fig. Al presents the emission of the three additional sources in our disk sample for the lines can be directly imaged: 
EC 82, R CrA and T CrA. LLN 19, a massive embedded object, also shows extended emission. As found previously 
for the embedded sources CSS 30 and IRS 43 by Herczeg et al. (2011), most of the visually extended sources appear 
to be outflows with the spatially extended emission occurring on only one side of the spectral trace, probably in the 
outflow cavity wall (see Figure Al). 

Evidence for outflows is also seen in broad absorption lines for these and other sources (Figures 3 and 4) (see also 
Thi et al. 2010; Herczeg et al. 2011). In contrast with the slow molecular winds analyzed through spectro-astrometry 
(Pontoppidan et al. 2011a), the outflows show large blue-shifts up to a few 100 km s~^ from source velocity. 
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Fig. A1. — Continuum subtracted 2-diniension stacked spectra of sources with spatially extended CO lines. The integrated flux spectrum 
is overplotted in white to show which components of the line profile are extended. EC 82 and T CrA are at P.A. 90, while R Cra and LLN 
19 are at P.A. 0. 



INNER RADII 

The gas closest to the star rotates most rapidly producing the wings of the line profile. Assuming that the CO lines 
arise from such a disk, the innermost extent of the CO can be measured based on the width of the line base. Two 
different methods of measuring the width of the stacked line profiles were used to overcome this problem (Table 8). 
In the first, the width was measured between the points where the line rose above 3 times the continuum noise. In 
the second, the width was measured at 5% of the peak flux values. These two measurements are often similar for our 
typical noise levels. When the continuum noise is low the first method provides a better grasp of the actual extent of 
the gas. Inner radii were determined by stacking clean lines to decrease continuum noise. 

The derived number is actually in terms of v sin i and thus the resulting radius is dependent on the inclination and 
stellar mass. The corresponding _Rsin i was calculated assuming Keplerian rotation and a typical stellar mass for the 
spectral type (O, B - 5 Mq, A - 3 M©, F, G - 1 M©, and K, M - 0.5 M©). The resulting numbers are highly dependent 
on inclination particularly for low inclination sources. For example, TW Hya, which has a known inclination of 4-7° 
(Pontoppidan et al. 2008; Qi et al. 2004), has an Rsin^i of 9.9 AU but an actual inner radius of less than 0.048-0.15 
AU. Also, any disk turbulence will broaden the line resulting in a smaller estimated inner radius than is actually the 
case. The disk turbulence is likely to be small compared to the line widths of >50 km s~^. 

The distribution of v sin i for the T Tauri stars is fiat within the errors as expected if inclination is the dominant 
effect. The range of values is consistent with Keplerian emission from the inner disk close to the dust sublimation 
radii as measured by infrared interferometry (e.g. Akeson et al. 2005). The smallest wsini values belong to the single 
peaked Keplerian sources, which in our sample are mainly transition disks with known low inclinations. The largest 
values tend to belong to the double peaked Keplerian sources. 

GALLERY OF SPECTRA 

Spectra of the entire CRIRES sample are presented in Figures 20-25. Figures 20 - 23 show 4.64-4.63 ^m [^^CO 
R(0)-P(7)] for all the disks in alphabetical order (Fig. 20: A-EL, Fig. 21: EX-IR, Fig. 21: IR-SY, Fig. 21: SY-Z). 
Eight disk sources, mainly in Chamaeleon, were observed only in the longer wavelength settings. These sources are 
denoted with * in Figs. 20-23 and the spectra are shown in Figure 24 covering 4.77-4.85 ^ [^^CO P(12)-P(18)]. The 
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22 embedded sources, including higher mass sources not covered in Herczeg et al. (2011), are shown in Fig. 25. 
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TABLE 1 
Summary of CRIRES disk observations 



Source 


Wavelength^ 
(nm) 


Date of observation^ 




AA Tau 


4616,4730,4868 


2007 Oct 12 




AS 205 N (A) 


4662.1,4676.1,4760.8,4773.6; 4730 


2007 Apr 21; 2007 Sep 1, 2008 May 2, 2007 Aug 
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AS 205 S 


4662.1,4676.1,4760.8,4773.6 


2007 Apr 21 




AS 209 


4716,4730 


2008 May 1 




CRBR 2422.8-3423 


4716,4730;4868 


2008 Aug 3; 2008 Aug 4 




CVCha 


4716,4800,4946.2 


2009 Jan 3 




CWTau 


4716,4868,4946.2;4730 


2008 Dec 30; 2009 Jan 1 




DFTau 


4716,4730 


2009 Jan 1 




DG CrA 


4716,4730,4833;4946.2 


2007 Apr 25; 2008 Aug 9 




DG Tau 


4716,4730,4868 


2007 Oct 15 




DoAr 21 


4716,4730 


2008 Apr 30 




DoAr 24E S (B) 


4716,4730;4730 


2007 Sep 2; 2007 Sep 3 




DoAr 24E N (A) 


4716,4730 


2007 Sep 2 




DoAr 44 (Haro 1-16) 


4716,4730;4833 


2007 Sep 4,2008 Apr 27;2008 Apr 29 




DRTau 


4716,4730,4833;4716,4946.2 


2007 Oct 10;2008 Dec 30 




EC 82 


4716,4730 


2007 Apr 22 




Elias 23 


4716,4730,4868;4946 


2008 Aug 06; 2010 Mar 4 




EX Lup 


4716,4730;4716,4730,4868 


2008 Apr 27; 2008 Aug 6 




FN Tau 


4868;4730 


2007 Oct 16; 2008 Dec 29 




FZ Tau 


4716 


2009 Jan 1 




GQ Lup 


4662.1,4676.1,4760.8,4773.6;4730 


2007 Apr 21; 2008 May 2, 2008 Aug 4 




HBC 680 


4716,4730 


2007 Sep 1 




HD 135344 B 


4716,4730,4929.3,4946.2 


2007 Apr 22 




HD 142527 


4716,4730,4868;4710;4076,4101 


2008 Aug 5; 2008 Aug 7; Aug 9 




HD 144432 S & N* 


4710,4730 


2008 Aug 2 




HD 144965 A (N) & B* 


4716,4730 


2008 Apr 29 




HD 176386 


4730 


2007 Sep 1 




Haro 1-4 A (S) & B* 


4716,4730,4833 


2008 Apr 29 




IM Lup 


4716,4730,4833 


2008 Apr 26 




IQ Tau 


4716,4946.2 


2008 Dec 30 




IRS 46 


4730 


2008 Aug 7 




IRS 48 


4730;4730,4833;4710,4730;4868 


2007 Sep 5; 2008 May3; Aug 2; Aug 5 




IRS 51 


4716,4730 


2008 Aug 3 




LkHa 330 


4730;4868;4716,4730 


2007 Oct 10; Oct 12; 2008 Dec 29 




RCrA 


4760.8;4730;4946.2 


2007 Apr 21; Sep 1; 2008 Aug 9 




RNO 90 


4716,4730,4833,4929. 3;4730 


2007 Apr 25; Apr 26 




RR Tau 


4730,4868 


2007 Oct 15 




RU Lup 


4716,4730,4833,4929. 3;4730 


2007 Apr 26; 2008 Apr 27 




RW Aur 


4868 


2007 Oct 15 




RYLup 


4716,4730,4840,4929. 3;4716,4730,4833 


2007 Apr 24; 2008 Apr 26 




S CrA N & S 


4716,4730;4730;4946.2 


2007 Apr 22; 2007 Sep 3; 2008 Aug 9 




Serp 32 


4716,4730 


2007 Sep 2 




Serp 64 


4716 


2007 Aug 30 




SO 411 


4868 


2007 Oct 14, Oct 16 




SR 9 A (S) & B* 


4716 


2008 Aug 8 




SR21 


4716,4730;4730,4833 


2007 Aug 30; Aug 31 




SR 24 A (S) & B/C(N) 


4716,4730,4833 


2008 May 3 




SX Cha A & B* 


4800,4946.2 


2008 Dec 30 




SYCha 


4800,4946.2 


2009 Jan 3 




Sz 68 A,B & C 


4716,4730,4833,4929.3 


2007 Apr 26 




T Cha 


4800,4820;4946.2 


2008 Dec 29; Dec 30 




T CrA 


4716,4730;4730,4833,4929.3 


2007 Apr 23; Apr 26 




T Tau N & S (A) 


4730 


2007 Oct 10 




TW Cha 


4800,4820,4946.2 


2009 Jan 2 




TWHya 


4662.1,4760.8,4773.6;4730 


2007 Apr 21; Apr 26 




TYCrA 


4716,4730 


2007 Apr 24 




UXTau 


4730,4868 


2007 Oct 12 




VSSG 1 


4716,4730,4868 


2008 Aug 6 




VV CrA N & S 


4716,4730,4840;4770,4779.5;4946.2 


2007 Apr 24; Aug 31; 2008 Aug 9 




VVSer 


4662.1,4676.1,4760.8,4773.6;4730 


2007 Apr 21; 2008 May 1 




VW Cha A & B/C* 


4800,4820;4946.2 


2008 Dec 29; Dec 30 




VZ Cha 


4800,4820;4946.2;4716 


2008 Dec 29; Dec 30; 2009 Jan 3 




Wa Oph 6 


4716,4730,4833 


2008 May 1 




WL 22 


4716,4730 


2008 Aug 3 




WX Cha 


4946.2;4830 


2008 Dec 30; 2009 Jan 2 





''This table summarizes the disk sources. See appendix table 6 for the complete embedded source list (also Herczcg 
et al. 2011 for low mass embedded sources). 

^Semicolons (;) divide different groups of observations (e.g. for VW Cha A, the 4800 and 4820 settings were both 
observed on Dec 29 while the 4946.2 setting was observed on Dec 30). 

^Serp 32 is a Spitzer selected source at 18:28:45.6 -00:07:21.6, and is probably a background giant/AGB star based on 
the high radial velocity. The source is not counted in the disk statistics. 

^Serp 64 is a Spitzer selected source at 18:29:01.7 -1-00:29:38.7. 
Detected companion was too faint for further analysis 
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TABLE 2 

Categorization of line profiles of disk sources 





1. Narrow 


2. Broad single peaked 


3. Broad 


4. Double 


No lines 


FN Tau 


As 205 N 


AS 209 


AATau 


HD 144965 A 


HD 135344 B 


DRTau 


CV Cha 


GQ Lup 


HD 176386 


LkHa 330 


S CrA N 


DF Tau 


IRS 48 


IMLup 


SR21 


S CrA S 


DC CrA 


RNO90 


SO 411 


TWHya 


RULup 


DoAr 24 E S 


VV Ser 


SY Cha 




VV CrA S 


DoAr 44 


WaOph 6 


UXTau 




VW Cha A 


EX Lup" 


WX Cha 






VZ Cha 


FZTau 
HD 142527 
HD 144432 S 
Haro 1-4 A 
IRS 51 
RRTau 
RW Aur 
SR9 A 
SX Cha A 
TWCha 
VSSG 1 






5. Narrow absorption 


6. Emission with broad 


7. Blue 


8. Blue emission/ 


Unclassifiable 


only 


central absorption 


absorption 


red absorption 




CRBR 2422.8-3423 


CWTau 


DC Tau 


AS 205 S 


Serp 32 


DoAr 21 


Elias 23*' 


HBC 680 


EC 82 


T TauN 


DoAr 24 E N 


IQ Tau 


IRS 46 


RCrA 




T Cha 


RYLup 


Sz 68 A/B 


SR24 A 




TYCra 


Scrp 64 


Sz 68 C 


SR 24 B/C 




WL 22 


TTau A 


VV CrA N 


T CrA 





''EX Lup was undergoing outburst during this time which affected the line profiles. A detailed 
examination of the CO emission based on this data can be found in Goto et al. (2011). 

''The CO v=2-l line profiles are double peaked so the underlying v=l-0 emission line profile could 
fall in category 4. 
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TABLE 3 
Radial velocities 



Source 


y*(Hdio) 


y*(LSR) 


Vco(LSR) 


Ref. 




(km s-i) 


(km s-i) 


(km s-i) 




AA Tau 


16.1 ± 2 


5.6 


9.0 


1 


AS 205 N 


-9.4 ± 1.5 


2.0 


5.0 


2 


AS 209 


-8.5 ± 1 


5.0 


5.7 


2 


CV Cha 


16.1 ± 1.3 


5.5 


3.0 


3 


CWTau 


14.5 ± 2 


5.6 


-8.7 


1 


DF Tau 


12 ± 5 


2.1 


8.7 


1 


DG CrA 


-1.8 ± 0.4 


5.4 


0.7 


3 


DGTau 


19.3 ± 2.7 


9.5 


-3.4 


4 


DoAr 24 E N 


-8.0 ± 2.5 


2.3 


3.5 


2 


DoAr 44 


-5.9 ± 0.4 


4.7 


11.4 


2 


EXLup 


-1.5 ± 2.3 


3.6 


3.3 


2 


FN Tau 


14.9 ± 0.4 


6.1 


4.6 


5 


GQ Lup 


-3.2 ± 0.6 


2.8 


6.0 


2 


Haro 1-4 A 


-7.6 ± 0.1 


2.8 


5.4 


2 


HD 135344 B 


-3 ± 3 


1.3 


7.4 


6 


RU Lup 


-0.9 ± 1.1 


4.6 


2.0 


2 


RW Aur 


16.0 ± 1.7 


6.5 


-50 


4 


RYLup 


-2.0 ± 0.8 


2.9 


-2.0 


2 


SCrA 


-0.3 ± 1.7 


7.1 


3.3 


2 


Sz68 


-2.3 ± 0.6 


3.9 


4.5 


2 


T Cha 


14.0 ± 1.3 


4.1 


6.0 


3 


TTauN 


19.1 ± 2 


7.6 


8.7 


1 


TWCha 


17.8 ± 1. 


7.0 


5.0 


3 


TWHya 


12.4 ± 0.2 


3.0 


3.3 


7 


VW Cha A 


17.2 ± 0.5 


6.6 


6.0 


3 


VZ Cha 


18.0 ± 1.1 


7.3 


12.7 


2 



^Herbig & Bell (1988); ^Melo (2003); ^Guenther et al. 
(2007); * White & Hillenbrand (2004); ^ MuzeroUe et al. 
(2003); ^ Dunkin ct al. (1997); ^ Setiawan ct al. (2008) 



TABLE 4 
Summary of ^^CO rotation diagrams 



Source 



T(i8CO) 
(K) 



Opt. thin 
radius (AU)* 



Nco 100AU2 
(10" cm-2) 



Gas mass 

(10-4 Mq) 



AS 205 N 


480 


0.8 


1.2 


0.0058 


DoAr 24E S 


490 


0.35 


0.089 


0.00045 


DRTau 


570 


1.3 


0.95 


0.0048 


HD 135344 B 


250 


5.0 


41. 


0.21 


HD 142527 


420 


3.1 


0.26 


0.088 


IRS 48 


270 


5.0 


40. 


0.20 


LkHa 330 


240 


13 


460. 


2.3 


RU Lup 


310 


4.4 


25. 


0.12 


S CrA N 


380 


5.6 


17. 


0.086 


S CrA S 


280 


8.9 


120. 


0.60 


SR21 


220 


18 


770. 


3.9 


TTau N 


640 


0.8 


0.54 


0.0027 


Where Acmit 


= -«?.nit 
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TABLE 5 

Vibrational flux ratios 



Source 


Vibrational fiux ratio 




(CO v=2-l/v=l-0) 


AA Tau 


0.1 


AS 205 N 


0.2 


AS 209 


<0.2 


CVCha 


0.25 


CWTau 


0.26 


DF Tau 


0.37 


DoAr 24E S 


0.12 


DoAr 44 


<0.068 


DRTau 


0.20 


FZTau 


<0.15 


GQ Lup 


0.19 


HD 135344 B 


0.067 


HD 144432 S 


<0.12 


IRS 48 


0.5 


LkHa 330 


0.11 


RNO 90 


0.18 


RU Lup 


0.29 


S CrAN 


0.42 


S CrAS 


0.2 


TWHya 


0.053 


WaOph 6 


0.25 


VSSG 1 


0.27 


VV CrA S 


0.4 


VVSer 


<0.2 


VW Cha A 


0.13 


VZ Cha 


0.31 



TABLE 6 

CRIRES OBSERVATIONS OF EMBEDDED YOUNG STARS 



Source 


Wavelength 




(nm) 


CrA IRS 2 


4716,4730,4833,4929.3 


EC 90 N & S 


4716,4730,4929.3,4946.2 


Elias 29 


4716,4730,4868 


Elias 32 


4716,4730,4833;4716;4730 


CSS 30 


4716,4730;5115 


HH 100 


4710,4730,4868,4946.2;4716,4730 


IRS 43 


4716,4730,4868;4831 


IRS 44 E & W 


4716;4730,4833;4946.2;4868;4716,4730 


IRS 63 


4716,4730,4833,4929.3;4716,4730,4868 


L1551 


4716,4730 


L1689 IRS 5 


4716,4730;4716,4730,4800,5115 


LLN8 


4716,4730 


LLN 17 


4716,4730;4716,4730,4868 


LLN 19 


4716,4730,4833 


RE 50 


4716,4730,4868;4770 


RNO 91 


4716,4730,4929.3,4946. 2;4770,4779. 5; 




4716,4730;4800,4820 


TMC lA 


4716,4730;4800,4946.2 


WL 12 


4716,4730;4833 


WL 6 


4716,4730,4833 



Date of observation-' 



2007 Apr 25 

2007 Apr 23 

2008 Aug 8 

2008 May 2; 2008 Aug 8; 9 

2010 Mar 4; 2008 Aug 4 

2008 Aug 2; 2007 Sep 2 

2008 Aug5; 2008 Aug 6 

2008 Apr 26; 29; 30; 2008 Aug 5; 6 

2007 Apr 25; 2008 Aug 5 

2007 Oct 16 

2010 Mar 10; 2010 Mar 15 

2007 Oct 11 

2007 Oct 11; 2007 Oct 13 

2007 Apr 25 

2007 Oct 11,2007 Oct 17 

2007 Apr 23; Apr 24; 

2010 Mar 19; 2010 Mar 9 

2010 Feb 10; 2010 Feb 3 

2007 Sep 1; 2007 Sep 3 

2008 Apr 26 



'^Equivalent to Table 1 for the embedded sample 
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TABLE 7 
Summary of observed line morphology^ 



Source 


M flux 

(Jy) 


P(8) flux 
(Jy kms-i) 


12CO 


i^CO 


i^CO 2-1 


AA Tau 


0.3" 


11.0 


E/A 


E 


E 


AS 205 N 


4.3'' 


112 


E 


E 


E 


AS 205 S 


5.4*' 


16.8/-59.7 


E/A 


A 




AS 209 


1.4= 


15.2 


E/FA 






CRBR 2422.8-3423 


0.4= 


-4.6 


A 


A 




CVCha 


1.4= 


32.9 


E/FA 


E/FA 


E 


CWTau 


1.6'' 


44.9/-14.3 


E/A 


E/A 


E 


DF Tau 


Lb-i 


43.2 


E/FA 


E 


E 


DG CrA 


0.5= 


7.0 


E/FA 






DGTau 


24d,e 


25.8/-74.3 


E/A 


E/A 


E 


DoAr 21 


0.88-^ 


-0.26 


A 


A 




DoAr 24E N 


0.15-^'= 


-0.04 


A 


A 




DoAr 24E S 


0.07^'= 


1.7/-0.1 


E/FA 


E/FA 


E 


DoAr 44 


0.54/ 


12.2 


E 






DRTau 


1.39 


43.8 


E 


E 


E 


Elias 23 


1.4= 


15.1/-22.3 


E/A 


E/A 


E 


EC 82 


0.6= 


17* 


E/A 


E/A 


E 


EX Lup 


1.3'' 


20.8 


E 


E 


E 


FN Tau 


0.26' 


1.3 


E 






FZTau 


1.0" 


51.1 


E 




E 


GQ Lup 


1.0= 


26.6 


E 




E 


Haro 1-4 S 


0.6= 


10.0 


E/FA 




E 


HBC 680 


0.5= 


-4.3 


A 


A 




HD 135344 B 


2.3^ 


6.4 


E 


E 


E 


HD 142527 


6.8^ 


40.9 


E 


E 


E 


HD 144432 S 


1.9= 


10.1 


E 


E 




HD 144965 A 


0.5= 


<0.73 


FA 






HD 176386 


0.4*= 


<0.33 


FA 






IM Lup 


0.6= 


<0.93 


FA 






IQ Tau 


0.5" 


3.3/-2.9 


E/A 




E 


IRS 46 


0.6= 


-16.1 


A 


A 




IRS 48 


2.0^ 


4.2 


E/FA 


E/FA 


E 


IRS 51 


0.92-^ 


12.1 


E/FA 


FA 




LkHa 330 


0.9' 


4.1 


E/FA 


E 


E 


RCrA 


69"^ 


-183/159/-375 


E/A 


E/A 




RNO 90 


1.3= 


44.7 


E/FA 


E/FA 


E 


RR Tau 


1.1™ 


1.7 


E/FA 






RU Lup 


1.1" 


32.1 


E 


E 


E 


RW Aur 


0.9= 


- 


E 






RYLup 


1.2= 


8.6/-3.0 


E/A 


E/A 




S CrAN 


1.1" 


42.5 


E/FA 


E 


E 


S CrAS 


3.3" 


74.3 


E/FA 


E 


E 


Serp 32 


2.8^ 


53.2/-70.4 


E/A 






Serp 64 


0.3/ 


4.0/-3.0 


E/A 


A 


E 


SO 411 


0.4= 


- 








SR9 A 


0.5/." 


5.3 


E 






SR21 


1.2° 


2.8 


E/FA 


E/FA 


E 


SR 24 A 


1.4/,c 


42.3/-28.9 


E/A 


A 


E 


SR 24 B/C 


0.6-'''= 


8.3/-6.9* 


E/A 


A 




SX Cha A 


0.6= 


8.9 


E 






SYCha 


0.5= 


<3.83 








Sz 68 A/B 


1.3P'= 


5.8/-6.3 


E/A 


A 




Sz 68 C 


0.1P'= 


8.3/-6.9* 


A 


A 




T Cha 


1.29 


-* 


A 


A 




T CrA 


3.7™ 


7.0/-18.7 


E/A 


A 




TTauN 


lid 


_* 


E 


E 


E 


TTauS 


15^* 


_. 


E/A 


E/A 


E 


TW Cha 


0.5= 


15.3 


E 






TWHya 


0.3= 


2.7 


E 




E 


TYCrA 


0.81™ 


-0.3 


A 


A 




UXTau 


3.9= 


_* 








VSSG 1 


i.o-f 


30.5/-1.2 


E/FA 


E/FA 


E 


VV CrA N 


07fe,c 


-7.8 


A 


A 


E 


VV CrA S 


2.2*='= 


26 


E/FA 


E/FA 


E 


VVSer 


3.0™ 


15.2 


E/FA 


FA 




VW Cha A 


0.4'' 


21.1 


E 




E 


VZ Cha 


0.8'" 


35.5 


E 




E 


Wa Oph 6 


0.9'' 


13.1 


E/FA 




E 


WL 22 


0.6= 


-1.0 


A 


A 




WX Cha 


0.6= 


- 


E 




E 



^E - emission, A - absorption, FA- foreground absorption, only present for 
J <8. Positive fluxes are in emission while negative fluxes denote absorption. 
Strongly affected by telluric features 
^a.) Luhman et al. (2006); b.) McCabe et al. (2006): c.) Derived from 
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TABLE 8 
Line widths 



Source 


width 3x noise 


width 5% peak 


Ri^sin'^i 


Inchn- 




(km s-i) 


(km s-l) 


(AU) 


ation(°) 


AA Tau 


172 


172 


0.060 


75 


AS 205 N 


103 


78 


0.29 


25 


AS 209 


165 


90 


0.22 




CVCha 


151 


154 


0.15 




CWTau 


142 


142 


0.088 




DF Tau 


102 


118 


0.13 




DG CrA 


200 


200 


0.089 




DGTau 


142 


155 


0.073 




DoAr 24E S 


139 


135 


0.097 




DoAr 44 


90 


102 


0.17 




DRTau 


62 


54 


0.62 


37 


EC 82 


74 


76 


0.30 




EX Lup apr 


173 


250 


0.028 




EX Lup aug 


134 


150 


0.079 




FZTau 


50 


105 


0.16 




GQ Lup 


139 


139 


0.091 




Haro 1-4 A 


64 


68 


0.38 




HD 135344 B 


27 


32 


3.43 


10 


HD 142527 


35 


40 


2.20 




HD 144432 S 


51 


64 


2.57 


48 


IRS 48 


17 


17 


37.4 




IRS 51 


66 


109 


0.15 




LkHa 330 


34 


35 


2.92 


10 


RNO 90 


186 


181 


0.11 




RULup 


135 


122 


0.12 




RYLup 


51 


146 


0.083 


85 


S CrAN 


101 


99 


0.36 




S CrAS 


166 


133 


0.10 




Scrp 64 


99 


123 


0.11 




SR9 A 


44 


141 


0.090 




SR21 


13 


16 


13.7 


20 


SR24 A 


122 


145 


0.085 




SR 24 B/C 


54 


83 


0.26 




T CrA 


39 


71 


0.70 




TTauN 


83 


146 


0.17 




TTauS 


62 


95 


0.39 




TWCha 


67 


132 


0.10 




TWHya 


12 


13 


9.9 


4 


VSSG 1 


118 


117 


1.31 


53 


VV CrA S 


60 


102 


0.17 




VVSer 


74 


79 


2.84 




VW Cha A 


98 


139 


0.091 




VZ Cha 


70 


99 


0.18 




Wa Oph 6 


114 


197 


0.046 


39 


WXCha 


127 


332 


0.016 





